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I.  INTRODUCTION 


The  Barrofcnnan  Method  [I]  provides  an  excellent,  low-cost,  tool  for  the 
evaluation  of  the  aerodynamic  model  of  a  rocket.  The  primary  limitation  is 
its  restriction  to  noncontrollable  sheet  fins.  This  report  presents  the  basic 
Barrowman  Method  which  pertains  to  the  body  of  the  rocket.  Also  presented  are 
modifications  to  the  basic  method,  which  allows  for  the  incorporation  of  air¬ 
foil-shaped  fins  with  control  surfaces  and  wing/fin  combinations. 

The  resulting  Modified  Barrowman  Method  (MBM)  was  applied  to  the  aft-wing 
configuration  of  PDAMS.  This  was  done  so  that  a  comparison  of  MBM  results  and 
wind  tunnel  data  could  be  made.  A  second  comparison  was  made  with  MBM  results 
and  the  current  PDAMS  aeromodel.  The  comparisons  are  presented  herein  for 
evaluation  of  the  method. 

An  aerodynamic  model  of  the  OPDAMS  vehicle  has  been  generated  through  the 
use  of  Che  MBM. 

II.  THE  MODIFIED  BARROWMAN  METHOD 

The  Barrowman  Method  is  a  means  of  calculating  the  aerodynamic  charac¬ 
teristics  of  a  rocket.  The  method  is  based  on  USAF  DATCOM,  with  simplifica¬ 
tions  of  fixed  sheet  fins.  The  method  is,  based  on  normal  force  coefficients 
and  is  only  valid  in  the  linear  regime.  The  modifications  to  the  method  allow 
for  the  use  of  published  data  for  airfoiled  fins  with  control  surfaces  and 
wing/fin  combinations.  The  results  of  interdigitation  of  the  wings  and  fins 
in  a  cruciform  configuration  is  demonstrated. 

A.  The  Barrowman  Method 

The  first  element  of  the  Barrowman  Method  is  that  the  normal  force 
contribution  of  a  straight,  constant  diameter  body  tube  is  zero.  Only  the 
nose,  body  diameter  transition  sections,  and  fins  contribute  to  the  normal 
force  of  the  rocket.  Calculations  are  performed  with  the  normal  force 
coefficients.  All  centers  of  pressure  (CP)  are  referenced  to  datum  zero  (0), 
which  is  located  at  the  nose. 

The  centers  of  pressure  of  nose  cones  and  body  transition  sections 
are  located  by  dividing  the  volume  of  the  component  by  the  area  at  the  com¬ 
ponent's  highest  datum.  This  gives  the  center  of  pressure  location  from  the 
component's  base.  A  transform  is  required  to  calculate  the  datum  of  the  CP. 

All  nose  cones  have  the  same  normal  force  coefficient  values,  provid¬ 
ing  the  base  area  is  the  reference  area.  Equation  (1)  gives  the  normal  force 
coefficient  for  nose  cones. 

(Cn„)n  -  2.0  (1) 

Convention  will  be  set  such  that  the  area  of  the  base  of  the  nose  is  the 
reference  area.  Otherwise,  a  scaling  factor  equal  to  the  nose  base  area 
divided  by  the  reference  area  is  required  to  multiply  (C^^  )n. 
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The  normal  force  coefficients  of  body  transition  sections  are  calcu¬ 
lated  by  Equation  (2), 


(C„^)t  -  2.0 


^2  V 


k'^ref , 


v'^ref 


(2) 


where  r^^  is  the  radius  of  the  body  at  the  end  of  the  transition  section  with 
the  lower  datum,  and  r^  is  the  body  radius  at  the  higher  datum  of  the  tran¬ 
sition  section.  The  reference  area,  S^gf,  is  equal  to  the  area  of  a  circle 
with  radius  r^gf.  It  can  be  seen  that  an  expansion  section  yields  a  positive 
(Cn^)T>  while  a  contracting  section  provides  a  negative 

For  a  continued  description  of  the  Barrowman  Method,  ic  will  be 
assumed  that  the  normal  force  coefficient  and  the  longitudinal  and  radial 
locations  of  the  center  of  pressure  of  the  fins  are  known.  The  means  of 
calculating  these  parameters  will  be  presented  in  modifications  to  the 
Barrowman  Method. 


The  total  normal  force  coefficient  of 
matlon  of  all  of  the  constitutant  normal  force 


the  rocket  is  obtained  by  sum- 
coefficients 


Cn^  ~  (3) 

The  location  of  the  overall  center  of  pressure  of  the  vehicle,  Z,  is 
given  by  Equation  (4). 

Z  -  Q(Cn^)iZi)/(r(Cn^)i)  (^) 


The  datum  of  the  center  of  pressure  of  the  normal  force  contributing  com¬ 
ponents  are  the  Z^. 

With  the  preceeding  data  the  moment  coefficients  can  be  calculated. 
The  reference  length  is  the  body  diameter  in  the  reference  section,  D  =*  2rjgf. 
The  0^2  coefficient  is  given  by  Equation  (5). 


(SM)  =■  (Z-W)/D  (6) 


The  vehicle's  center  of  gravity  is  given  by  W.  The  static  margin  is  denoted 
by  (SM).  Thus,  for  a  stable  r£c^et  design  the  static  margin  must  be  greater 
than  zero,  or  in  other  words,  Z>W. 

The  damping  moment  coefficient,  is  given  by  Equation  (7). 

Cm^  “  I  (Zi-W)2/D2  (7) 
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Since  the  nose  and  body  transition  sections  exhibit  radial  symmetry, 
they  do  not  contribute  to  the  roll  moments.  Thus  the  fins  are  the  only  con¬ 
tributor  to  the  rolling  moments.  The  radial  location  of  the  CP  of  the  fins, 

Y  ,  is  found  by  calculation  of  the  radial  centroid  of  the  exposed  fin  area, 
r 

It  is  assumed  that  the  vehicle's  center  of  gravity  lies  on  the  vehicle's 
centerline.  The  calculation  of  the  C^p»  the  rolling  moment  coefficient  due  to 
a  roll  Mte,  P,  is  calculated  by  Equation  (7),  where  Y.  is  substituted  for 
Zi  and  W  is  the  radial  location  of  the  center  of  gravityi  multiplied  by  the 
number  of  wing  pairs  (two,  for  cruciform  wings). 

The  constants  are  used  in  Equation  (8-10),  which  are  the  homogeneous 
equations  of  motion  for  a  non-thrusting  rocket. 

I.  a  +  ^  u  Sref  D  Cjj*  a  +  ^  u  D  a  “  0  (8) 

I^S  +  I  u  Sref  6  +  ^  u2  S^gf  D  -  0  (9) 

P  +  P  u  Sref  P  “  0  (10) 

Due  to  the  symmetry  of  the  rocket,  the  a“derlvatives  and  8“derivatives  are 
equal. 

For  vehicles  that  are  propelled  by  rockets,  the  exhaust  jet  of  the 
engine  provides  a  damping  term  that  needs  to  be  included  in  Equation  (8)(9). 

The  Cnj.  is  given  by  Equation  (11)  and  its  use  is  shown  in  Equation  (12).  For 
completeness  of  the  method,  similar  augmentation  of  the  g  equation  is  required. 

C®**"  «  i(Ln-  -  w)^/D^  (11) 

“oi 

exh 

1^0  +  £.  u  Sygf  Cth.o  +  C  a  +  Sj-gf  D  o  “  0  (12) 

L  2  ®  2  ® 

Where  m  is  the  mass  flow  rate  of  the  rocket  engine  and  L^g  is  the  datum  of  the 
rocket  nozzle's  exhaust  plane.  At  engine  burnout,  m  becomes  zero  and  Equation 
(12)  becomes  Equation  (8). 

B.  Modifications  of  the  Barrowman  Method  for  the  Use  of  Airfoiled  Fins 
with  Control  Surfaces  and  Airfoiled  Wings 

The  techniques  developed  for  the  modifications  to  the  Barrowman  Method 
have  been  generalized  and  can  be  applied  to  a  wide  variety  of  configurations  and 
airfoils. 
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To  start,  obtain  the  two-dimensional  airfoil  data  from  a  source  such  as 
Hoerner  [2]  or  Riegels  [3].  Since  rockets  generally  use  symmetric  airfoils  for 
fins,  »  0  and  Cm  =  0.  For  a  non-symmetric  airfoil  the  non-zero  Cn  and 

u  o  j 

Cjj  terms  can  be  added  to  the  following  and  terms.  The  transformation  of 
the  infinite  span  lift  coefficient,  (C|^  )o,  to  a  finite  span,  Cj  ,  is  given  by 
Equation  (13).  ^ 

(Czjw  -  [1  +  (Cj^)o  (13) 

Where  e  *  0.95  for  elliptic  wing  loading.  The  lift  coefficient  is  now  area 
scaled  to  the  reference  area  of  the  rocket.  Considering  a  wing  with  constant 
chord,  c,  and  span,  S,  the  wing  area,  A^,  is  the  product  of  the  chord  and  span. 
The  virtual  airfoil  passing  through  the  rocket's  body  is  required  for  these 
calculations.  Usually  there  are  some  cut-outs  in  the  airfoil  for  control  boxes 
and  such  that  must  be  taken  into  account.  This  is  accomplished  by  subtracting 
the  cut-out  area,  A^,  from  the  wing  area  A^.  Equation  (14)  presents  the  wing 
lift  scaled  to  the  rocket's  reference  area. 

*  (^w  ~  ^c)(C£^)w/Sref  (14) 

With  the  lift  acting  at  the  quarter  chord,  and  since  these  calculations  are 
valid  in  the  linear  regime,  the  small  angle  approximation  can  be  applied.  This 
suggests  that  the  normal  force  coefficients  are  equal  to  the  lift  coefficients. 
At  this  point,  the  Barrowman  Method  calculations  for  static  stability  coeffi¬ 
cient  can  be  done. 

The  control  surface  calculations  are  also  started  from  a  source  such  as 
Hoerner  or  Riegels  [2,3].  These  sources  contain  wind  tunnel  data  for  various 
control  surfaces  on  different  airfoils.  The  data  indirectly  gives  Cjl5  and 
as  the  data  presented  is  usually  9a/96  and  For  Cjjj  the 

chain  rule  is  applied  to  the  result  of  Equation  (14). 


(Cjg)fs 


Is. 

35 


(15) 


The  control  surfaces  are  usually  not  full-span  due  to  such  things  as  control 
boxes  and  not  extending  to  the  wing  tips.  To  account  for  this,  an  empirical 
method  based  on  Weisslnger's  Method  [4]  is  used.  The  graphs  allow  for  the 
determination  of  Kf  and  Kf^.  The  factors  are  calculated  from  these  parameters 
determined  from  graphs  [4].  The  result  of  Kf  or  Kf^^  is  dependent  on  whether  the 
outer  or  inner  span  of  the  control  surface  is  used.  If  the  control  surface 
extends  to  the  wing  tips,  then  Kf  *  1. 


Kf(Kfj^) 


-  k^^l  +  k^  (/?• 
+  kj  sin^tan"^ 


(tan 


(16) 


*  (*^f  ■  J^f f )(Cjig)fs 


(17) 
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For  5^0  there  is  a  moment  about  the  quarter  chord.  Equation  (18) 
is  a  means  of  calculating  the  torque  couple  about  the  wing's  quarter  chord. 


■'m 


(c/4)6 


(Kf  -  Kf .) 


c  ^w~^c  ^^m(c/4) 
^  ^ref  3'S 


(18) 


This  term  is  added  to  the  torque  term  generated  by  the  lift  of  the  control 
surface . 


(w  r 

D  C 


^6  ^m(c/4)6 


(19) 


This  generates  a  forcing  term  for  Equation  (8)  as  shown  in  Equation  (20). 

5e  (20) 


=  !.  Sref  D 
2  ■ 


Since  and  0^  are  both  negative,  this  implies  that  positive  pro- 

^  (c/4)5 

duces  a  negative  pitching  (yaw)  moment.  A  sign  change  may  be  performed  at 
Cnjg  to  agree  with  the  convention  used. 


The  roll  derivative  C,  is  calculated  from  C/..  and  the  radial  cen- 

Lda 

troid  of  the  control  surface  (Yj,)cs*  Since  each  control  surface  is  indepen¬ 
dent  of  the  other,  the  aileron  control  surface  angle  (fi^)  is  obtained  from 
unequal  elevator  or  rudder  flap  pair  deflections. 


written 


oe  *  (^2  +  d4)/2  (21) 

6r  -  (6i  +  53)/2  (22) 

5a  -  (62  -  54  +  5^  -  (23) 

With  the  calculation  of  ,  the  nonhoraogeneous  roll  equation  can  be 

^«a 


(Yp)cs  C  /O 
®a 

^  u  Si-gf  D  C  6a 
2  ^5a 


(24) 

(25) 


This  completes  the  Modified  Barrowman  Method  for  a  standard  rocket 
configuration  such  as  PDAMS.  For  the  configuration  of  OPDAMS,  (basically  a 
PDAMS  with  forward  wings)  the  effects  of  the  wings  on  the  fins  need  to  be 
known.  The  primary  effect  is  downwash  induced  by  the  wings  on  the  tail. 
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For  wings  Chat  are  in  line  with  the  fins,  the  calculation  of  the 
downwash  is  straightforward  [4].  For  interdigitated  wing/fins  the  calcula¬ 
tions  need  to  be  modified.  Due  to  symmetry,  consider  Che  wings  to  be  level 
with  some  angle  of  attack  to  Che  airflow.  For  cruciform  wings  and  fins,  Che 
plane  of  Che  fins  is  inclined  at  7r/4  rad  for  interdigitation.  Thus,  projec¬ 
tion  of  Che  fin  areas  onto  the  plane  of  the  wing  is  given  by  Equation  (2b). 

Ap  =  2A.pC0S  ^  (26) 

The  angle  of  attack  generates  two  components  from  the  rotation  of  the 
fins  to  an  interdigitated  position.  One  component  is  spanwlse  and  the  other 
is  in  the  plane  perpendicular  to  the  span.  The  projection  of  the  airflow  onto 
the  plane  perpendicular  to  the  flow  yields  an  effective  angle  of  attack  for 
the  fins  of 

E 

Op  “  Ctcos  Tr/4  (27) 

E 

The  lift  of  the  fins  becomes  Equation  (28): 


T  1  2 

L  =■  —  0  a  C»  u  A_c(. 
2  a  ?  . 

*  T  Cj  Aa 


(28) 


For  calculations.  Equation  (28)  proves  that  cruciform  fins  inter¬ 
digitated  with  cruciform  wing  by  it/4  rad  can  be  considered  Co  be  in  line,  with 

9e » 


one  exception  -  downwash.  The  downwash,  is  performed  in  the  same  manner, 

3a  — . 

but  now  the  aerodynamic  center  is  above  the  plane  of  the  wings  by\  2  Vj./2, 
This  vertical  displacement  reduces  the  effect  of  downwash.  The  downwash  is 
reduced  by  31  percent  for  the  OPDAMS  vehicle.  Thus,  the  fins  experience  a 
greater  angle  of  attack,  which  Increases  lift  and  eliminates  the  requirement 
of  larger  fins.  Therefore,  the  advantage  to  interdigitation  is  significant. 
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III.  PDAMS  Aerodynamic  Models 


Two  PDAMS  aerodynamic  models  have  been  formulated  via  MBM.  The  first 
model  corresponds  to  the  aft-wing  configuration  used  by  Killough  during  wind 
tunnel  testing  [5].  A  comparison  of  the  MBM  theoretic  model  and  the  wind  tun¬ 
nel  data  is  presented  as  a  demonstration  of  the  technique.  An  aeromodel  of 
the  current  PDAMS  vehicle  was  calculated  from  the  elements  of  the  aft-wing 
configuration.  This  was  compared  to  the  current  PDAMS  aeromodel.  There  was 
good  agreement  of  the  MBM  model  and  the  wind  tunnel  data.  There  were  dif¬ 
ferences  with  the  current  aeromodel,  with  the  major  difference  being  in  Che 
Cnj^  term.  The  implications  of  this  major  difference  are  discussed. 

A.  Comparison  of  the  MBM  Model  with  Wind  Tunnel  Data 

The  PDAMS  aft-wing  configuration  aerodynamic  coefficients  were 
derived  from  wind  tunnel  data.  The  data  was  presented  in  tabular  form  for  a 
range  of  a  and  6*  The  derivatives  were  calculated  by  averaging  the  slopes  of 
the  data  in  Che  linear  regime,  -0.175  rad  ^  o  ^  0.175  rad.  This  was  done  to 
correct  a  6  percent  asymmetry  in  the  data  about  a  *  0.  Due  to  the  symmetry  of 
Che  vehicle,  full  data  for  5  was  not  taken,  but  rather  only  for  5  0  rad. 

Thus,  the  delta  derivatives  were  calculated  for  only  fi  £  0.  The  wind  tunnel 
data  is  presented  by  the  averages,  with  the  variances,  both  positive  and 
negative,  shown  in  parentheses. 


TABLE  1.  Comparison  of  MBM  Model  with  Wind  Tunnel  Data 


Wind 

Tunnel 

MBM 

Percent 

Coefficient 

Data 

Results 

Difference 

18.96 

+0.94 

-1.12 

19.04 

0.4 

8.10 

+0.63 

-0.26 

8.15 

0.6 

16.92 

+3.14 

-2.23 

19.88 

15 

Cmj 

-15.14 

+3.66 

-0.44 

-14.11 

7 

8.05 

+0.51 

-0.43 

8.82 

9 

Z(in) 

27.65 

+1.40 

-0.92 

28.56 
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B.  Comparison  of  the  Current  PDAMS  Aeromodel  with  MBM  Results 

The  current  PDAMS  configuration  is  not  the  same  as  the  vehicle  that 
was  tested  in  the  wind  tunnel.  Thus,  as  a  second  demonstration  of  MBM,  it  was 
applied  to  the  current  configuration.  Configuration  data  was  obtained  from 
MICOM  Drawing  RLC-2077.  The  current  aeromodel  was  obtained  from  Mr.  Sam 
Stauffer  [6]  of  Rexham  Aerospace. 

TABLE  2.  Comparison  of  Current  PDAMS  Aeromodel  with  MBM 


Coefficient 

Current 

Aeromodel 

MBM 

Results 

18.03 

19.04 

7.84 

8.15 

9.13 

8.57 

-11.87 

-9.29 

52.03 

3a/35 

1 

• 

O 

-1.08 

Static 

Margin 

0.51 

0.45 

The  normal  force  coefficients  and  roll  moment  coefficients  are  the 
same  for  both  the  current  model  and  the  aft-wing  wind  tunnel  model  for  MBM. 
This  is  in  agreement  with  the  results  of  Killough  [5]  for  the  aft-wing  and 
forward-wing  wind  tunnel  configurations.  But  there  is  a  difference  in  the 
normal  force  coefficients  between  the  current  configuration  and  the  wind  tun¬ 
nel  model. 


The  most  significant  difference  with  the  current  PDAMS  aeromodel  was 

in  Cm  .  To  demonstrate  the  importance  of  this  term,  a  pull-up  maneuver  was 

0 

formulated.  This  was  done  by  setting  6  ■  0.262  rad  (15  )  with  the  vehicle 
traveling  in  the  standard  sea-level  atmosphere  with  a  constant  velocity  of 
150  m/ s  (approximately  80  percent  of  terminal  velocity.)  The  initial  calcula¬ 
tions  were  performed  in  a  zero  gravitational  field.  The  current  PDAMS  aero¬ 
model  has  a  turning  radius  of  370  m  while  the  MBM  model  has  a  turning  radius 
of  502  m.  The  quarter  turn  where  the  vehicle  enters  on  the  downward  vertical 
and  the  impact  is  at  the  bottom  of  the  loop.  A  gravitation  field  was  superim¬ 
posed  on  the  loop,  stretching  it  out.  The  current  aeromodel  was  used  to  set 
the  entry  and  Impact  point.  The  MBM  aeromodel  was  put  through  the  same 
maneuver  with  the  same  entry  point.  The  result  was  a  miss  of  the  impact 
point  by  168  m. 

This  quarter  loop  maneuver  is  an  analytic  example  to  indicate  the 
effects  of  a  reduced  C^j  for  a  system  with  a  heavy  emphasis  on  the  end-game 
guidance  strategy. 


8 


IV.  OPDAMS  Aerodynamic  Model 


Two  changes  to  the  Baseline  OPDAMS  vehicle  were  made  to  accommodate 
information  received  during  analysis.  First,  from  the  PDAMS  calculations,  it 
was  observed  that  there  was  a  difference  in  the  span  of  the  fins  between  PDAMS 
and  OPDAMS  Baseline.  Since  the  OPDAMS  aft-end  is  a  PDAMS  aft-end,  in  order  to 
eliminate  the  need  to  change  the  fin  spans  on  the  aft-end,  the  standard  PDAMS 
fins  were  placed  on  OPDAMS.  This  also  allows  use  of  PDAMS  wind  tunnel  data  to 
finalize  the  aerodynamic  model. 

The  second  change  was  due  to  a  shift  of  the  CG  by  1.05  inches  rearward  of 
the  design  location  [7].  This  produced  a  3a/36  *  6.4,  which  is  unacceptable. 
An  equation  was  formulated  which  locates  the  CG  as  a  function  of  the  quarter 
chord  of  wing  datum  ,  Z^. 

W  =  23.369  in  +  0.1035  Z„  (29) 

This  was  based  on  a  full-up  OPDAMS  vehicle  with  aluminum  wing  ring  [8].  The 
equations  of  Cm  (4)  and  Cmr  (19)  were  written  as  functions  of  Z^.  This 

Qt  6 

system  of  equations  were  solved  simultaneously  for  3a/36  “  1.7.  The  result 
was  Zy,  =  20.5,  the  leading  edge  of  the  wings  at  19.0.  The  result  of  these  two 
changes  were  designated  OPDAMS  Baseline  Block  II. 

TABLE  3.  OPDAMS  Baseline  Block  II 


Length  -  51.25  in 

Diameter  -  6.00  in 

Wing  Span  -  27.00  in 

Wing  Chord  -  6.00  in 

Wing  Leading  Edge  -  19.00  datum 

Fin  Span  -  20.00  in 

Fin  Chord  -  9.00  in 

Fin  Leading  Edge  -  42.25  datum 

Surface  -  33.3  percent  chord,  overhanging  nose  balanced  flap 
pivoted  about  flap's  quarter  chord,  each  flap  set  independently. 

Airfoils  -  NACA  0012,  cruciform  arrangement  with  constant 
chord  and  no  sweep,  wings  and  fins  interdigitated. 

Mass  -  60.89  lb 

Ir  -  725  lb  in2 
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TABLE  4.  OPDAMS  Baseline  Block  II  Aerodynamic  Model 


36.234 

8.15 

‘^ma 

15.959 

-27.027 

183.096 

^LP 

110.482 

8.829 

3a/36 

-1.694 

W 

25.491  datum 

Z 

28.1337  datum 

Static  Margin 

0.44 

^long 

133,000  lb  in2 

C  1.18  (analytic  estimate) 


V.  CONCLUSIONS 


The  Modified  Barrowman  Method  has  been  presented  and  compared  to  wind  tun¬ 
nel  tests.  The  agreement  of  MBM  results  and  wind  tunnel  data  was  excellent. 

The  MBM  was  applied  to  OPDAMS  and  the  aerodynamic  model  for  OPDAMS  Baseline 
Block  II  was  presented. 

While  the  theoretic  results  of  MBM  characterized  the  wind  tunnel  results, 
it  must  be  stressed  that  they  do  not  replace  wind  tunnel  data.  The  MBM  is  an 
excellent  tool  for  analysis  of  rocket  designs.  There  are  pathologic  cases  for 
which  any  analytic  method,  including  MBM,  will  fail.  It  is  the  wind  tunnel, 
an  Immense  analog  computer,  which  can  best  indicate  the  true  aerodynamic  model. 
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